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We study the behavior of a Bose-Einstein condensate in which atoms are weakly coupled to a 
highly excited Rydberg state. Since the latter have very strong van der Waals interactions, this 
coupling induces effective, nonlocal interactions between the dressed groundstate atoms, which, 
opposed to dipolar interactions, are isotropically repulsive. Yet, one finds partial attraction in 
momentum space, giving rise to a roton-maxon excitation spectrum and a transition to a super solid 
state in three-dimensional condensates. A detailed analysis of decoherence and loss mechanisms, 
suggests that these phenomena are observable with current experimental capabilities. 
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Since introduced by Landau in a series of seminal arti- 
cles [l] , the notion of a roton minimum in the dispersion 
of a quantum liquid has been pivotal to understanding 
superfluidity in helium. This later led to the prediction 
of a peculiar solid state upon softening of the roton ex- 
citation energy [2], simultaneously possessing crystalline 
and superfluid properties. In such a supersolid the 
particles that must supply the rigidity to form a crys- 
tal, at the same time provide for superfluid nonviscous 
flow. Forty years after its conjecture, this apparent con- 
tradiction continues to attract theoretical interest and 
has ushered in an intense search for experimental evi- 
dence in solid ^He, whose interpretations are currently 
under active debate [4, 5]. 

Here, we demonstrate how three-dimensional roton ex- 
citations can be realized in atomic Bose-Einstein con- 
densates (BECs), thereby introducing an alternative sys- 
tem to study supersolidity. The supersolid phase transi- 
tion is shown to arise from effective interactions, realized 
through off-resonant optical coupling [6-8] to highly ex- 
cited Rydberg states. Owing to the strong increase of 
atomic interactions with their principal quantum num- 
ber n, resonantly excited Rydberg gases have proved 
to be an ideal platform to study strong interactions in 
many-body systems [qI on short /is time scales. The 
present approach — based on off-resonant two-photon 
excitation [see Fig(T^] of Bose-condensed alkaline atoms 
— permits to utilize the strong Rydberg interactions over 
much longer times of ~ 100ms. In particular, we consider 
coupling to nS Rydberg states with vanishing orbital an- 
gular momentum, which, as opposed to dipole-dipole in- 
teractions, gives rise to isotropically repulsive interaction 
potentials for the groundstate atoms, and, thus, ensures 
stability of the condensate. 

The system is described as a gas of N atoms with 
mass M at positions r^, each possessing a ground state 
\gi) and an excited nS Rydberg state, denoted by \ei). 
The two states are optically coupled with a two-photon 
Rabi frequency Vt and detuning A (see Fig(T^). Defin- 
ing corresponding transition and projection operators 
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FIG. 1: (a) Schematics of the considered three- level atom, il- 
lustrating the laser coupling between the atomic groundstate 
\noS) and the Rydberg state \nS). For Ai ^ Qi, the sys- 
tem reduces to an effective two-level atom, with the states 
1^) = \noS) and |e) = \nS) coupled with a two-photon Rabi 
frequency Q, and detuning A. (b) Effective potential result- 
ing from the off-resonant coupling to the strongly interacting 
Rydberg states for n = 60 and A = 50MHz. Panels (c) 
and (d) provide an enlarged view of the potential showing 
the contributions from both groundstate-Rydberg atom and 
groundstate-groundstate atom interactions (solid line) as well 
as the sole contribution from the latter (dashed line). 



(<^7 /5 = e,g), the resulting A/'-particle in- 



teraction can be written as 



i<j 



where Hl = ^ Xli ^eg + o-ge describes the laser cou- 
pling and Vee{Tij) = Ce/rfj > denotes the van der 
Waals (vdW) interaction between two Rydberg atoms at 
a distance Vij = — vj. Because of the strong Cq ~ n^^ 
scaling of the vdW coefficient, such Rydberg- Rydberg 
atom interactions are orders of magnitude larger than 
those of groundstate atoms. We are interested in the 
potential surface 1^0(1*1, tat) that asymptotically con- 
nects to the many-body groundstate \G) = \gk), un- 
der the condition of far-off resonant driving ^^/|A| <C 1 
and for A < 0. Under the latter condition the A^-body 
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FIG. 2: Dispersion relation e{k) for different values of the 
interaction parameter a. The arrow indicates the roton gap 
S. 



potential energy of the many-body groundstateIG) is 
separated from the excited states by at least A [10|, jus- 
tifying the use of a Born-Oppenheimer (BO) treatment 
for the atomic dynamics on the potential surface Wg- 
The BO surface Wg can be determined from a many- 
body perturbation expansion in the small parameter 
r^/|A| <C 1, which up to fourth order only involves cou- 
pling to the singly and doubly excited many-body states 



= \ei)(S)k^i \9k) and \Eij) 
plicitly, we obtain 



eiej)(S)k^ij \9k)- Ex- 



Wg 



E 



i 



\{G\Hi^\E,)\^ 
HA 



AHA 



{2nA-Vee{r^J)) 



IQTi^A^ 



Kj 



(1) 



Omitting constant terms, the leading-order, A/'-body BO 
surface is given by a sum of binary effective potentials 
^{^ij) — C'e/l^fj + ^c)' exemplarily shown in Fig(Tl3. 
The pecuHar shape of U is easily understood within a 
simple two- atom picture: For far-distant atomic pairs 
(vij ^ i?c), a small fraction (r2/2A)^ is indepen- 
dently admixed to each groundstate atom, such that 
U{rij) is of vdW-type with an effective coefficient Cq = 
(r^/2A)^C6. At smaller distances, however, the inter- 
action shift V^e(rij) renders dressing to the doubly ex- 
cited \Eij) states ineffective, such that the effective po- 
tential approaches a constant value below a critical dis- 
tance Rc = (C6/2n|A|)i/^ Typically, can take on 
rather large values of a few /im. 

For simplicity, the above discussion has ignored inter- 
actions between pairs of groundstate atoms as well as 
groundstate-Rydberg atom interactions. Both, however, 
can be included in the calculations of the effective po- 
tential (see Figdt). Although being comparably strong, 
the range of the corresponding contributions is consider- 
ably smaller than Rc. Upon avoiding photo-excitation of 
Rydberg- molecular resonances [11] these additional po- 
tentials may, thus, be described in terms of an s-wave 
scattering pseudopotential, with a scattering length a. 
On the other hand, the large value of Rc prevents such 



a simplified treatment for U{r). However, as long as 
Cq < n^Rt/M its effect can be described within a first 
Born approximation 12 1 [26|. With these simplifications 
and in the zero-temperature limit one arrives at the fol- 
lowing nonlocal, nonlinear Schrodinger equation 
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for the condensate wavefunction ^(r). The contact in- 
teraction term (x g = AnTi^a/M may be tuned to very 
small values. For clarity, it will, hence, be omitted in the 
following and briefly discussed at the end of this work. 

For small Cq the BEC ground state corresponds to a 
homogeneous superfluid with density po- Its elementary 
excitations with wave number k and corresponding en- 
ergy e are calculated from the corresponding Bogoliubov- 
de Gennes equations, which yield e(/c)^ = \^[\^ + 
2poU{k)], where U{k) is the Fourier transform of the in- 
teraction potential U{r). Upon appropriate length and 
energy scaling (see Fig 12]), the dispersion is determined 
by a single dimensionless parameter a = pqMCq/Ti?Rc^ 
parametrizing the interaction strength. Asymptotically 
the dispersion relation e(k) has phonon and free particle 
character, at small and large /c, respectively. However, 
due to the inner potential plateau at r < Rq^ the momen- 
tum space potential U has negative attractive contribu- 
tions around k — krot ^ 27r/i?c, such that the spectrum 
develops a roton minimum ai k = krot (see Fig|2|) for suf- 
ficiently large a > Ofrot ^4.8. The corresponding roton 
gap 6 decreases with increasing a and ultimately vanishes 
at Qfinst ^ 50.1 (see FigllK), marking the onset of a roton 
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FIG. 3: (a) Roton gap ^ as a function of the interaction pa- 
rameter a. (b) Energy density e = |^(^| ~fj'' + ^|^) 
for different crystal symmetries relative to the energy den- 
sity Ehom = Ti'^a/S of a homogeneous BEC. Panel (c) pro- 
vides an enlarged view around the transition point. The 
effective Rydberg state lifetime for excitation of ^^Rb to 
n = 60(C6 = 9.7 • 10^°a.u. 21] ) with A = 50MHz and 
po = 10"^^cm~^ is shown in (d). 
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FIG. 4: Snapshots of the BEG dynamics for a time- varying 
interaction parameter a{t) shown in (a). Panels (b)-(e) show 
the density along orthogonal slices through the simulation 
box at times indicated in (a). The upper and right axes in 
(a) show the actual time and Rabi frequency for a ^^Rb BEG 
with n = 60, A = 50MHz and po = 2 • lO^^cm"^. 



instability, at which density modulations may develop 
without energy cost. In ultracold e: ases, a similar behav- 
ior was found in tightly confined, pancake-shaped dipolar 
BECs [13, J^]. While it is the partially attractive nature 
of the anisotropic interactions that generates the 2D ro- 
ton excitations in this case [15], it also leads to collapse in 
the post-roton-instability phase, thus, precluding forma- 
tion of supersolid density modulations 16|, [l7| in higher 
than one dimension [18]. The present interactions, on 
the other hand, are entirely repulsive, which may indeed 
provide for stable regular density modulations. 

Addressing this question, requires us to go beyond the 
Hnear analysis discussed above. We performed numerical 
simulations of Eq.(|2]) on a three-dimensional grid of ~ 10^ 
grid points with periodic boundary conditions. Because 
of the nonlinear character of the Gross Pitaevskii equa- 
tion a standard imaginary-time integration may generally 
not converge to the true groundstate of the system. As 
also found for lattice-confined dipolar bosons [19], there 
is a large number of stationary states, that correspond to 
local minima of the total energy. As a consequence, the 
system generally approaches a glassy state with short- 
range ordered density modulations, when starting from 
a homogenous initial state (see below). We, thus, used 
variational calculations, based on periodically arranged 
Gaussians with varying width and lattice constant, to 
provide the proper initial wavefunction for a subsequent 
imaginary time-evolution according to Eq.(|2]). 



Some of the obtained energies are shown in Figl3]3. For 
small values of a the BEG groundstate is a homogenous 
superfiuid. At a critical value of aguso ^ 30.1, one finds 
a transition to a stable supersolid state. This first-order 
transition precedes the roton-instability [20] and takes 



place at a finite roton gap of Ss 
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istence of several competing states with similar energies 
but different crystal symmetries (see FiglSb), may gener- 
ally complicate the experimental preparation of ordered 
states. In this respect, the dynamical tunability of the 
interaction strength Cq via changing the laser intensity 
can serve a useful tool to steer the BEG evolution. 

In order to demonstrate this point, we also studied 
the time evolution, starting from a homogenous BEG. 
As a specific example, we discuss the BEG dynamics for 
a simple time-dependence of a, shown in FigH^. The 
calculation starts from a homogenous condensate with 
small random phase noise and uses a complex-time in- 
tegration with a small imaginary contribution [2^. The 
instantaneous increase of a at time t = from a = 
to 60 induces the roton instability. This sudden param- 
eter quench, however, causes relaxation towards a short- 
range ordered, "glassy" [5] state (FigHjD), as discussed 
above. As a is decreased close to the phase transition 
some of the structures vanish entirely, leading to a mixed 
phase in which extended superfiuid fractions of nearly 
constant density coexist with density-modulated domains 
(Fig lit). The latter increase in time (FiglHi), and ul- 
timately merge to form sizable "crystallites" of regular 
density modulations (FigH^). 

Turning to a discussion of the experimental feasibility, 
we consider a particular example of coupling to 60S Ry- 
dberg states in a ^^Rb condensate, for which Rydberg 
excitation has recently been demonstrated [23]. Figl5] 
shows the corresponding "phase diagram" for a typical 
density of lO^^cm"^, also including a finite s-wave scat- 
tering length a > 0. The latter only leads to some in- 
crease of the critical Rabi frequencies for inducing the 
roton instability. Importantly, the transition to a su- 
persolid can be realized with Rabi frequencies of a few 
hundred kHz, and the condition |A| ^} can be well 
fulfilled deep in the roton-instability regime. Yet, the de- 
tuning is sufficiently small to avoid excitation of nearby 
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FIG. 5: Required laser parameters to induce the roton insta- 
bility (solid lines) for n = 60 and po = lO^^cm"^. The upper 
curve also accounts for s-wave scattering with a = lOOao- 
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Rydberg states and near-resonant dipole coupling to ad- 
jacent pair states for distances ^ Rc. For typical values 
of i?c7 the number of atoms within a single density peak is 
on the order of 10^, which justifies the applied meanfield 
description in terms of Eq.(|2]). 

Major limitations for the stability of Rydberg gases 
generally stem from the finite lifetime of the involved 
excited states and from auto-ionization of close Ryd- 
berg atom pairs, initiated by near-resonant dipole-dipole 
couplings to energetically close pair states [2^. In 
the present case, the latter are suppressed by a factor 

(2a) (r^ +i^g)^ interaction blockade of dou- 

bly excited states, and, thus, limited to negligible small 
values of SOmHz for the parameters range of Figl3]i. Ryd- 
berg state decay, due to spontaneous emission and black- 
body radiation, is more significant with typical rates on 
the order of 7 ~ 10ms~^. However, since the Rydberg 
state is only weakly populated by the far off-resonant 
coupling to the groundstate, the resulting effective decay 
rate 7eff = (r^/2A)^7 can be decreased to much smaller 
values. As shown in Figl3]i, for typical atomic densi- 
ties and laser parameters the effective lifetime is as large 
as several seconds, over the entire range of relevant in- 
teraction parameters. A second loss mechanism arises 
from spontaneous decay of the intermediate P state (with 
rate 7'), used to drive the 5S n'P nS two-photon 
transition (see Fig(T^). Again, the corresponding effec- 
tive decay rate 7^^ = (l^i/2Ai)^7' is greatly suppressed 
for small ratios l^i/Ai, which, however, also reduces the 
two-photon Rabi frequency Q = l]il]2/2Ai. With Q2 
of a few hundred MHz and Qi/Ai = 10~^ one realizes 
a two-photon Rabi frequency of several hundred kHz, as 
considered in the example of FiglSl Such large Rydberg 
Rabi frequencies (1^2), could be achieved by choosing a 
6P state for the two-photon transition, allowing to drive 
the Rydberg transition with a strong 780nm-laser. At 
the same time this gives a longer intermediate-state life- 
time 7'"^ = 121ns [25], and with Qi/Ai = 10"^ yields 
long effective lifetimes of 7^^^ = 0.5s. Together with 
the small value of 7eff this ensures sufficiently low loss 
rates to observe the phenomena discussed in this work 
(cf. FiglH). 

In summary, we have shown that off-resonant Ryd- 
berg excitation of ultracold atoms provides a promising 
route for the preparation of three-dimensional supersolid 
phases in Bose-Einstein condensates. The unambiguous 
realization of this elusive state opens up a range of new 
studies, from investigations of thermal effects to super- 
solid formation in finite systems and their dynamical re- 
sponse to perturbations, such as trap rotations. 

Generally, the described setting offers new avenues for 
the realization of complex nonlocal media, where the 
sign, shape and strength of interactions can be tuned 
through proper choice of the addressed Rydberg state 



and the applied laser parameters, and even permits local 
spatial control of nonlocal interactions via tightly focused 
beams. 
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